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ABSTRACT 

We derive the luminosity-temperature relation for the super-critically accreting black holes 
(BHs) and compare it to the data on ultraluminous X-ray sources (ULXs). At super-Eddington 
accretion rates, an outflow forms within the spherization radius. We construct the accretion 
disc model accounting for the advection and the outflow, and compute characteristic disc 
temperatures. The bolometric luminosity exceeds the Eddington luminosity i^Edd by a loga- 
rithmic factor 1 + 0.6 Inm (where m is the accretion rate in Eddington units) and the wind 
kinetic luminosity is close to LEdd- The apparent luminosity for the face-on observer is 2-7 
times higher because of geometrical beaming. Such an observer has a direct view of the inner 
hot accretion disc, which has a peak temperature T^ax of a few keV in stellar-mass BHs. The 
emitted spectrum extends as a power-law Fe oc E^^ down to the temperature at the spher- 
ization radius Tgp « m^^/^ keV. We associate Tmax with a few keV spectral components 
and Tsp with the soft, 0.1-0.2 keV components observed in ULXs. An edge-on observer sees 
only the soft emission from the extended envelope, with the photosphere radius exceeding the 
spherization radius by orders of magnitude. The dependence of the photosphere temperature 
on luminosity is consistent with that observed in the super-Eddington accreting BHs SS 433 
and V4641 Sgr. Strong outflows combined with the large intrinsic X-ray luminosity of the 
central BH explain naturally the presence of the photoionized nebulae around ULXs. An ex- 
cellent agreement between the model and the observational data strongly argues in favour of 
ULXs being super-critically accreting, stellar-mass BHs similar to SS 433, but viewed close 
to the symmetry axis. 

Key words: accretion, accretion discs - black holes physics - X-ray: binaries - X-ray: galax- 
ies 



1 INTRODUCTION 

A large number of ultraluminous X-ray sources (UL Xs) has been 
discovered in the nearby star-forming galaxies (see iMushotzkvl 
[2004, for a review). Their luminosities exceed considerably 2 x 
lO'^^erg s~^, the Eddington limit for a stellar-mass (~10Mq) 
black hole (BH). This was argued to be the evidence for ex- 
istence of th e intermediate-mass BHs ( IMBHs, with M — 
lO^-lO^M©; IColbert & Mushotzkvl [19991) . The large apparent 
luminosities can also be prod uced b y super-critical accretion 



on to a stellar-mass BH (jS hakura & S unvaevI 
SS73; Ijaroszvhski etai]|l98o[ lAbramowicz et al.l 



1973 , hereafter 
19881 : iLipunoval 



* E-mail: juri.poutanen@oulu.fi (JP); galja@sai.msu.ru (GL); fab- 
rika@sao.ru (SF) 



1 19991: lFabrikall2004 ISegelman et al.ll2006l). or by the geometric 
Jpabrika & Meschervakov' '200 1': 'King et al.l lioOll) or relativistic 
(Revnold s et al..,1997. : Kording et al. 2002j) collimation of radia- 
tion. 

A strong argument in favour of IMBHs is the presence of a 
soft, 0.1-0.2 keV comp onent in their spectra ( Kaaret et al. 200^ 
iMiller et all2003Ll2004l) . Arguments against the IMBH interpreta- 
tion i nclude theoretical problems with their formation l iKing et^ 
l200lh an d their non-standard spectra, which show a cutoff at a 
few keV (Stob bartetai] |2006'), while other BHs, stellar as well 
as super-massive, at a few pe r cent of Eddington l uminosity have 
hard power-law-like spectra jZdziarski et al .111 9971) . Furthermore, 
the observed anti-correlation between luminosity and temperature 
contradicts the L oc law expected for standard accretion discs 
l lFeng& Kaaret I2OO7I) and some ULXs show a harder, 1-4 keV 
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ther mal component with the corresponding radiu s of only 30-40 
km jMakishima et alj|2000l : IStobbart et alj|2006l) . All this raises 
further doubts on the IMBH interpretation. The soft components 
can be interpreted as signatures of an extended photosphere (SS73; 
lLipunovall99^ : lKing & Pounds 2003), but then we need to explain 
the simultaneous presence of the harder emission. 

The low-frequency, 0.02-0.2 H z, quasi-periodic oscilla- 
tions detected from ULXs in M82 (fstrohmaver & Mushotzkvl 
l2003h. Holmberg IX ([Pewangan et al ] |2006|) . and NGC 5408 
I Strohmaver et alj|2007[) . if interpreted as the Keplerian frequen- 
cies at the innermost stable orbit around a BH, argue in favour 
of IMBHs. However, oscillations with very similar frequencies 
ha ve been observed fro m the BHs in the Milky W ay, Cygnus X - 
1 jVikhlininetalll 19941) and GRS 1915-1-105 CMorgan et alj|l997h . 
which certainly are not IMBHs. 

Important clues on the nature of ULXs come from the pres- 
ence of extended photoionize d nebulae ar ound them ( Wang 2002; 
iPakuU & Mirionil2003i : lKaaret et alj2004h . The observations imply 
a rather isotropic source of the ionizing radiation rejecting the idea 
of a strong beaming of radiation. The nebulae are dynamically per- 
turbed with the velocity gradients of ^ 50-100 km on the scale 
of 50-100 pc (Leh mann et al. 2005 : iFabrika & Abolmas"ovll2006l : 
iRamsev et alJuOOfil) . This points towards the activity of the central 
engine in form of a wind or a jet. The ULX nebulae are similar to 
W50, the nebulae around SS 433, the only known persistent super- 
critical accretor in the Milky Way, radiating presumably around 
10'"^ erg in the UV iDolan et al. 1997). A super-critically ac- 
creting compact sourc e can produce strong winds and jets which 
can inflate the nebulae jBeg elman et al.*198d: 'Lehmann e t alj2003: 
Pakull et alj|2006l: Ipabrika & Abolmasov 2006; Abolma sov et al.l 



2007h . On the other hand, there is no good physical reason why 



a sub-critically accreting IMBH would produce such a strong out- 
flow. 

The observed similarities between W50 and ULX nebulae lead 
us to consider seriously the idea that the central engines of ULXs 
are super-critically accreting stellar mass BHs similar to SS 433. 
We do not see directly the X-ray source in SS 433, but if obs erved 
along the symmetry axis, it would be a bright X-ray source iKat 



at3 
Bvl 



1986), which we interpret as an ULX (Fabrika & Meschervaki 
,2001; King 2002; Fabrika 2004; Begelman et al. 200l). 

In this paper, we develop a model for the super-critical accre- 
tion disc accounting for the effects of advection and outflows. We 
also construct a one-dimensional, vertically integrated model of the 
wind, estimate the optical depth through the wind and determine 
characteristic temperatures as a function of the mass accretion rate 
and the luminosity. We further compare the resulting luminosity- 
temperature relations to the data on ULXs as well as the BHs in 
our galaxy and LMC. 



2 SUBCRITICAL ACCRETION DISCS 

The standard accretion disc theory (SS73) can be applied when the 
accretion rate is not very high and the luminosity does not exceed 
the Eddington limit 



^Edd = 



GAfA/Edd 
2i?i„ 



— 1.5 10 m 



1.7 



ergs ,(1) 



where A^Edd = A?,hGM/ck = 2 10^**m g s"^ is the Eddington 
accretion rate, 7?in = 3-Rs is the inner disc radius, Rs = 2GM/c? 
is the Schwarzschild radius. The stellar mass, measured in solar 
masses, is m = AZ/Mq, k = 0.2(1 + X) = 0.34 cm^ g~^ is the 



Thomson opacity and X is the hydrogen mass fraction (which we 
assume equal to solar X — 0.7). 

The energy flux from one face of the disc (in Newtonian ap- 
proximation) as a function of radius R is 



-(7?) 



3 GMM 
8^ R^ 



.-1/2 



(2) 



where Af is the accretion rate and r = R/Rin- This results in the 
effective temperature variation T{R) oc R~^^'^ at large R and the 
maximum observed color temperature 

Tcmax = 1.26/c m'^^^m^^ 



keV. 



(3) 



which is reached at rmax = (3/2)''/^ « 1.38. Here m = 
M / M-Edd is d imensionless accretion rate . The color correction 
/c ~ 1.7 (e.g. IShimura & Takahara|[T99^ describes the harden- 
ing of the spectrum relative to the black body. The emitted flux 
Fe oc at _E < fcTcmax. The total emitted luminosity. 



L = 



poo 

/ Q+{R)'iTiRAR = mLEdd, 



(4) 



depends on the maximum temperature as L oc T): 



i4 

c,max" 



3 SUPERCRITICAL ACCRETION DISCS 

The matter supply to the BH in a close binary can significantly 
exceed the Eddington rate, mo = Mo/M^dd ^ 1. There are two 
different views on the way the accretion proceeds in this regime. 

(i) The 'P o lish d oughnut' or slim disc models 
jjaroszvhski et al.l Il980l ; lAbramowicz et al] Il988l) assume 
that all the supplied gas reaches the BH, but most of the gravita- 
tional energy released in the disc is advected into the hole as the 
photons are trapped in the flow. Inside the trapping radius _Rtr 
(where photon diffusion and accretion time-scale are equal), the 
vertical component of gravity scales as R~^ (as the disc relative 
scale-height H/R is of the order unity) resulting in the same radial 
dependence of the radiation flux Qrad, because the disc is radiation 
pressure supported. This implies the effective temperature distribu- 
tion _R2ll^f_2Sii5SLLi'^ EFe) emitted spectrum Fe oc E^^ (see 
e.g. IWataraietalJIioog) . Naturally, a logarithmic dependence of 
bolometric luminosity of inner supercritical disc on Mo is yielded: 



fltr 



fitr 



Lx J Q,^d{R)RdR (X J i?-Mi? = ln^^^ «lnmo.(5) 



The last relation follows from the fact that the trapping radius scales 
with the accretion rate as Rti ~ moftn ^ Rm- The total disc 
luminosity exceeds the Eddington one by a logarithmic factor: 



ibol 



(1 + In mo). 



(6) 



(ii) Alternatively, instead of spending most of the dissipated 
energy to increase the entropy of the gas, this energy can be spent to 
eject the excess mass (SS73). In this model most of the gas is blown 
away by the radiation pressure, the accretion rate decreased linearly 
with radius M{R) oc R, and only a small fraction of it, m ~ 1, 
makes it to the holeQ The linear dependence of the accretion rate 



1 A .similar physics might operate at a very low accretion rate, when the gas 
cannot cool resulting in a positive Bernoulli pai'ameter. This implie s that the 
gas is efFectively unbound and can easily produce outflows i Naravan & Yil 
ll994lBlandford & Begelmarj|l999l) . 
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on radius results in the dependence of tiie radiative flux[f|The 
luminosity and the disc effective temperature radial distribution are 
very similar to the slim disc case (see below). The largest difference 
is in the presence of a strong wind which blocks and reprocesses the 
radiation from the central part of the disc. The observed spectrum, 
therefore, depends on the velocity structure and geometry of the 
wind and the position of the observer relative to the disc rotational 
axis. 

Reality might be somewhere in between, with both advection 
and o utflow operating together as shown by the numerical simula - 
tions lE ggum et alj|l988| - lohsuga et al.ll2005l : lokuda et alj|2005h . 
Below we first construct a one-dimensional, vertically integrated 
model of the supercritical disc with the outflow and then discuss 
the effect of the wind on the emergent spectrum. 

3.1 Supercritical disc with outflow and without advection 

The scale-height of accretion discs for mo ~ 1 is determined by 
the balance between the radiation pressure force and the vertical 
component of gravity (SS73), 

H . 3 



n . 6 \ 



1 - r 



-1/2 



(7) 



When the accretion rate exceeds the critical value rhcr = 9/4, the 
ratio H/R exceeds uni ty at some radii and the outflow is inev itable. 
Accurate calculations jSisnovatvi-Kogan & Blinnikovlll977h show 
that the outflow can start even at a smaller rate. At mo ^ mcr, the 
disc starts 'feeling' that it is supercritical at the spherization radius 
Tap = Rsp/Rin ~ rho, where the outflow starts. 

The wind affects the disc structure by removing the angular 
momentum. The modified angular momentum conservation equa- 
tion is (see e.g. Lipunova 1999.) 



dR\ ' ' J dR di? ' 

where luk{R) = {GM / R^Y^'^ is the Keplerian angular velocity. 



(8) 



g{R) = 2-K Tr4.R' 



(9) 



is the torque and Ti-^ ~ "iH tr4, is the vertically integrated viscous 
stress. The first term on the right hand side of equation ^ is re- 
sponsible for transfer of the angular momentum by viscosity, while 
the second term is the momentum transfer by the outflowing matter. 
The viscous heating rate per unit area is 



2Q+{R) = -Tr4,R 



dujvi 



3 giR) 



(10) 



dR 2 

If advection is not accounted for, all dissipated energy is con 
verted locally to the radiation: 



Q+{R)=Qr.i{R)- 



(11) 



We can assume that a fraction e„ of the radiative energy goes to 
accelerate the outflow: 



1 1 

ewQrad(^) ~ — 



dM{R) iuiR^ _ GM dM{R) 



2 2tvR dR 



87ri?2 dR 



(12) 



where v'^{R)/2 — lo^R? /2 i s the energy requ ired to eject a unit 
mass from radius R to infinity. ILipunoval ( 1 19991) assumed the max- 
imum outflow rate with ew = 1 and obtained an analytical solu- 



tion of the system of equations lf8l)-l ll2l (. applying boundary condi- 
tions M(i?sp) = Mo (where Rsp is a parameter of the model) and 
g{Rin) = 0: 



M(r) 

MEdd 



go 



"^0 T^VJ' ^^^=P' 



rsp 1 + ^r-p' 



(13) 



mo, 



r > rs 



mor 



3/2 



-5/2 



3rs 



go 



+ mo 



.1/2 



1/2 



r < rs. 



r > Ts 



(14) 



where go = MEddvGMftT. 

At intermediate radii 1 <C r < rgp, the accretion rate is close 



to the dependence A/(r) 
radiative energy flux is 



Mo r /rsp suggested by SS73 and the 



Qrad(i?) 



3 (GMy/^gir) 



GMM(R) „_2 
oc R . 



(15) 



8tt i?V2 8n R3 

Note that in the standard model there is a factor of 3 present in a 
similar formula (see equation |2ll, which is related to the fact that 
the energy dissipated close to the inner disc radius is transported to 
larger radii and radiated there. In the model with the mass loss, this 
factor is missing, implying that the energy transfer from the inner 
part of the disc is negligible. 

The luminosities produced within and outside rsp are 



L{r < rsp) 

L^dd 

L{r > rsp) 



mol'^''=P 



1 + fr.-//^ 



5 mo 



T- 1 I 2 -5/2 ■ 

■5 rsp -rs, 



3'sp 



Defining the spherization radius by the condition (SS73) 

L{r > rsp) = LEdd, 

we obtain 

5 . 

rsp « - mo, 

which gives the total luminosity 
3 



■ In mo 



(16) 



(17) 



(18) 



(19) 



and the Eddington accretion rate to the BH, m(7?in) = 1. A frac- 
tion of this luminosity can escape as radiation and the rest as kinetic 
energy of the outflow. 

It is worth noting here that if we assume in this model that 
only a small fraction of the dissipated energy is used to produce the 
outflow, i.e. Ew < 1 (see eq. I12t . then the resulting accretion rate 
is too high. With the high radiative efficiency (because we neglect 
advection) the locally emitted flux is strongly super-Eddington, and 
therefore the model is unphysical. However, accounting for the ad- 
vection allows us to construct a self-consistent model even in this 
case (see below). 



^ We note here that any supercritical disc model (with advection, with out- 
flows, or both) predicts such a dependence, which just follows from the 
vertical balance of gravity and radiative pressure. 



3.2 Supercritical disc with advection 

At supercritical accretion rates, the advective transport of the vis- 
cously generated heat to the black hole becomes important in the 
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1 10 lO"^ 10^ 10* 1 10 lO"^ 10^ 

Radius y—R/R^^ Radius r=R/R^^ 



Figure 1. (a) The radiative energy flux emitted from the disc (times r^, i.e. energy flux per logarithm of radius) as a function of radius for rho = 10 and 10"^. 
The flux is measured in units Qo = GMMEdd/STri??^ = c^/36GMk = 1.52 lO^'' m ^ erg cm s ^ . The solid curves from the top to the bottom are the 
exact solutions of the advective disc equations for ew = 0, 1/2 and 1. Bold dots represent the positions of the spherization radius defined by condition (TT). 
The dashed curves correspond to the analytical solution given by equations )10t . (TT) and )14t . (b) The accretion rate as a function of radius for mo = 10"^. 
The solid curves are the exact solutions for the advective disc, while the dashed curve is solution (13) for the non-advective disc with the outflow. The dotted 
curves are the linear approximations given by formula (25). 



energy balance equation. The diffusion time for photons, traveling 
to the disc surface, becomes larger than the characteristic time of 
the radial displacement of the matter. Due to the advective remov- 
ing of the heat, the locally radiated flux becomes smaller: 
1+ 



Qrad — Q Qadv ■ 



(20) 



A model of a supercritical advective disc with an out- 

flow was first proposed by iLipunovj il999 ) (se e also 
iKitabatake. Fukue & Matsumotol |2002| ; iFukud |2004) . The 
maximum mass loss from an advective disc by the energy-driven 
wind is about twice smaller than in the model without advection 
and amounts to about 3/ 5 of initial accre tion rate Mq. We modify 
the model considered by Lipunoval ( 1 19991) assuming that a fraction 
e„ < 1 of the radiation energy flux is spent on the production of 
the outflow (see eq. 1121 ). We solve the standard set of equations 
for the advective disc (as described in details in lLipunovij[T999l) 
for various mo and e„. 

The outflow occurs within the spherization radius rsp, which 
is defined self-consistently from condition ( 117b . For different tw, 
the spherization radius varies because a different amount of angu- 
lar momentum is gone with the wind, and, consequently, the struc- 
ture of the outer subcritical disc is different as it depends on the 
boundary condition at rsp. The results of our calculations can be 
approximated (with the accuracy of 2 per cent for rho > 5) by a 
simple formula: 



^ 1.34 - 0.4ew + O.le^ - (1.1 - 0.7e„ 

mo 



-2/3 



(21) 



The resulting radiation flux and the mass accretion rate as 
functions of radius are shown in Fig. [T] One sees that the radia- 
tive energy flux depends very weakly on Ew This also means that a 
similar solution will be obtained if a smaller mass is ejected with a 
larger velocity. The radial dependence of Qrad{r) is similar to that 
in discs with strong mass loss and no advection (see section |3T| and 
the dashed curve in Fig. [TJi). At r < rsp ~ mo, (5rad(r)r'^ is al- 
most constant. This is just the consequence of the fact that the disc 
is close to the Eddington limit locally everywhere, i.e. H/Ris close 
to unity. At r ^ rsp, Qrad decreases as r~^ as in the standard disc. 



The total luminosity is well approximated by equation l |19l l. As we 
assume that a fraction ew < 1 of the radiative energy flux is spent to 
drive the wind to infinity, a fraction 1 — Ew can escape as radiation. 
We note that in t he advective disc (with or without th e outflow) the 
scale-height (see lBeloborodovll 19981 : iLipunovall 19991 , and Fig.[TJi) 



H/R'. 



Qo 



< 0.8, 



(22) 



which justifies our use of the vertically integrated quantities. 

Because the outflow is optically thick at most radii and the 
radiation is partially trapped (see belo w), it is energy-driven (not 
momentum-driven as was assumed by iKing & Poundsll2003l) . and 
its kinetic luminosity can exceed I/Edd- 

The fraction of the initial accretion rate that passes through the 
inner radius can be approximated for rho > 2.5 as 



min 
rho 



1 - a 



1 



a rm-o 



-1/2 ' 



(23) 



where (0.83 — 0.25ew). For example, for mo = 1000 and 

Ew = 1/2 we get rsp = 1.16 rho and rh-in = 0.66 rho, while in the 
case of the maximal outflow with = 1 we have rsp = 1.04 rho 
and rhin = 0.43 rho- The total outflow rate for large rho is thus 



Mw ^ aM, 



0- 



(24) 



The exact solutions for the accretion rate inside the spheriza- 
tion radius can be approximated by the linear relation (compare 
solid and dotted curves in Fig.[Tj5) 



rh(r) ^ Thin -f (rho — hiin) — . 



(25) 



This can be easily understood from equation il2i : the radial deriva- 
tive of the accretion rate is proportional to r^Q{r), which is almost 
constant at r < rsp . The linear behaviour of the acc retion rate on 
radiu s is also obtained in the numerical simulations jOhsuga et al.l 
l2005h . 
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Figure 2. The Thomson optical depth through the wind as a function of 
radius in the direction parallel and peipendicular to the disc for niQ = 10^. 
The wind parameters are /3 = f = 1 and e^, = 1/2. The perpendicular 



optical depth reaches the maximum of ~37TI| 



1/2 , 



-^50 at the spherization 



radius rg 



' mo. 



4 STRUCTURE OF THE OUTFLOW 
4.1 Optical depth through the outflow 

The gas ejected from the accretion disc at cylindrical radius R 
gains the velocity perpendicular to the disc (in the z direction), 
Vz — (,vk{R) with C > 1. A mass element reaches asymptotically 
the velocity of (vk (R), with ( — y/^'^ — 1. If the angular momen- 
tum is conserved, the gas moves (in the ballistic approximation) 
along the line z/R = (at large radii. The radial (i.e. projected to 
the disc plane) velocity is thus (vk{R)/(, - Because the disc scale- 
height H/Rr^O.6 and ^ is not expected to exceed 1.5-2 (owing to 
the energy constraints), ("~0.6-1.7. The outflow is thus confined in 
the region outside the cone of opening angle 6 given by cot 9 — 
and it occupies r2„/47r = cos 6 =0.5-0.85 fraction of the sky. 

Let us construct the vertically averaged wind model. As our 
baseline we take the advective disc model with the outflow de- 
scribed in the previous section. We take the accretion rate as given 
by equation J25l > with rsp and min given by equations ( 121b and J23l >. 
The mass outflow rate (at r < rsp) is then 



Mw(i?) = 



dM{R) 
dR 



dR = AfEdd(rho — mi„ 



(26) 



Within the spherization radius the mean wind radial velocity should 
scale with the local Keplerian velocity. At i? > Rsp, the mass loss 
rate is constant and the gas reaches asymptotically velocity of ~ 
/3 vk{Rsp), where /3 ~ ~ L We can then approximate the wind 
radial velocity profile by a simple function: 




,(i?)= , 

P\/ -^77 , R > Rsp- 
From the mass conservation law, 
Mw(i?) ^C4ttR^ p{R) w„, 



(27) 



(28) 



we find the mean density p and the optical depth in the perpendic- 
ular direction, t±(R) = KpC,R: 



rj_(r) 



To mo — min 

13 rsp 



r < Ts 



, 1/2 -1 

(rsp - Ijrsp r , r > n 



(29) 



Here ro=MEdd\/6'v/47rci?in ~5. The optical depth has the maxi- 
mum rj_,max~3(my^-mo ^''^)e„//3 at r = rsp (see Fig.|2j. 

The Thomson optical depth from radius r in the direction par- 
allel to the disc is 



/■oo 1 r-oa ^ I 

r||(r)=«/ p{R')dR'^^ r^{r') — . 

J R S Jr ^ 



(30) 



The maximum of ~ 8(riT,|!,''^ — 4/3)ew/C/3 is reached at r = L At 
r > rsp it decays in the way identical to t±_ (Fig.|2}: 

ry (r) = rx(r)/C oc r"'^. (31) 

The outflow becomes optically thick at mo ~ 2.5 for ew ^ 1/2. 

4.2 Photospheres and the emitted spectrum 

We can define three characteristic radii (see Fig. [2}: (1) the radius 
of the inner photosphere rph,in. where rx = 1; (2) the spherization 
radius rsp ~ rh.o, where the optical depth through the wind in the 
normal direction is maximal; and (3) the outer photosphere rph, 
where the wind becomes transparent rx ~ T|| = L 

The inner photosphere is almost independent of rho: 

rph,in~l + /-. (32) 
The outer photosphere (for mo 2> 1), 



rpii 



o . 3/2 

3^mo , 



(33) 



is much larger than the spherization radius. 

A face-on observer would see the emission from three separate 
zones defined by the three characteristic radii: 



r < rph.in, 

rph. in ^ r < rsp, 

rsp < r < rph, 



zone A, 
zone B, 
zone C. 



(34) 



The characteristic disc temperatures can be obtained from the 
Stefan-Boltzmann law 



Qrad(ii) =<TSBr^(i?). 

For the advective disc with 
perature is about 



(35) 

1/2 the maximum effective tem- 



1.6 m 



^1/4 



1 - 0.2m, 



-1/3 



keV. 



(36) 



It is reached at rmax ~ 1-06 < rph, in for large mo and varies 
little with the accretion rate. Variations in tw affects this temper- 
ature only by a per cent or so (see Fig. [T^). The observed color 
temperature differs from Tmax by a color correction factor /c. The 
exact values for Tmax and r^iax, however, depend on general rela- 
tivity corrections, which are neglected here. In zone A, the wind is 
transparent (i.e. it is momentum-driven) and the radiation escapes 
unaffected by the outflow. 

In zone B, the wind is opaque and the energy generated in 
the disc is advected by the wind. The ratio of the photon diffusion 
time in the wind, t±(R/c, to the dynamical time-scale, (R/v^,, is 
Sew/v^ ~ 1 (this supports the view that the wind here is energy- 
driven). Thus the radiation escapes at a radius about twice the en- 
ergy generation radius. This does not change the radial dependence 
of the effective temperature T oc R^^^^, resulting in a power-law 
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Figure 3. (a) The luminosity-temperature relation for sub-critically accreting BHs in the Milky Way and Large Magellanic Cloud (data from lGierlinski & Pond 
[2004). Most objects (except LMC X-1) show a well pronounced correlation L oc T^-^^^ consistent with the standard accretion theor y (SS73). The bolometri c 
luminosity and the color temperature Tc.max are coiTected here for the effect of inclination and relativistic effects (see details in Gierliriski & Donell2004l) . 
Theoretical dependences are shown for 1.5, 5 and 10 solar mass BHs. (b) The luminosity-temperature relation for super-critically accreting BHs. The curves 
are theoretical dependences for the model of the advective disc with the outflow given by equations 1191 and )36M38t . The accretion rate is mo ~ 2 at the 
point of separation of various temperatures, and the curves continue until mo = 10"^ (filled circles at the curves indicate a change in mo by a factor of 10). 
The following temperatures are shown (from the right curve to the left): the maximal color disc temperature Tc.max = /cTmax with the color correction 
factor /c = 1.7, the color temperature at the spherization radius Tc,ap (also with /c = 1.7), and the temperature at the outer photosphere Tp^^. The wind 
parameters are /3 = = 1, Ew = 1/2. The upper (solid) set of curves is for the BH mass m = 15 and the lo wer dashed curves are fo r m = 5. The stars 
show the positions of the super-critically accreting stellar-mass BHs, SS 433 <DolanetalJll997h and V4641 Sgr <Revnivtsev e7d]|2002h . Other symbols are 
the apparent bolometric lumin osities and the tempe ratures obtained from the spectral fits with the blackbody and DISKBB model to the XMM-Newton data of 
a set of ULXs (from Table 5 of lStobbart et alj2006l) . The crosses show the data for LMC X-I if one increases the apparent luminosity by a factor of 5. 



spectrum Fe oc E ^ extending from about Tph,in to the tempera- 
ture at the spherization radius rsp ~ mo : 

Tsp ^ 1.5 m-^/^mo ^''^ (l + 0.3 niQ ^'"') keV. (37) 

The resulting temperature should also be reduced by (1 — tw)^''*, 
because some energy is transferred to the outflow. 

The outer zone C emits about the Eddington luminosity which 
is produced mostly in the disc at radii r > rsp. The photon diffu- 
sion time here is smaller than the dynamical time, thus most of the 
radiation escapes not far from the radius it is produced. This results 
in the effective temperature variation close to r""^/* and the nearly 
standard spectrum Fe oc E'^^^ (SS73). 

An edge-on observer would see only the blackbody-like emis- 
sion corresponding to the temperature at the outer photosphere, 
which for mo 3> 1 (i.e. rph 3> ^sp) takes the form: 

Tph « 0.8 (^^^ m-i/^rho keV. (38) 

For accretion rates slightly exceeding the Eddington, rph < rsp, 
and the dependence of Tph on mo is much stronger. 

At intermediate inclinations, the central hot part of the disc 
may be partially blocked by the wind, and an observer would see a 
soft spectrum peaking at Tsp. 



5 COMPARISON WITH OBSERVATIONS 

The standard model for sub-critically accreting BHs (SS73; Sect. 
O predicts the relation L oc T^^^^ oc M. At super-Eddington ac- 
cretion rates, three characteristic temperatures are identified: (i) the 
maximal color disc temperature Tc,max = /cTmax ~ 1.6/cm^^'^'' 
keV, (ii) the color temperature at the spherization radius Tc,sp ~ 
1.5fcm~^^*mQ^^^ keV, and (iii) the outer photosphere temper- 
ature given by equation l l38t . The bolometric luminosity and the 
temperatures depend parametrically on mo according to equations 
l |19l l and l|36|l-l|38t. These theoretical dependences are shown in 
Fig. nib). The luminosity observed along the symmetry axis may 
exceed iboi by a factor 1/(1 — cos 6')~2-10 for the outflow height 
z / R^Q.6-2. On the other hand, a fraction ew is spent on accel- 
eration of the outflow. Together these effects result in the 10-30- 
fold excess over L-^dd at high accretion rates (se e also lOhsuga et al.l 
l2005l : lBegelman et alJlOOalFabrika et al.l2006l) . Thus the absolute 
maximum apparent luminosity in our model is about 10^^ erg s"'^ 
for a 20 Mq BH. 

The BHs in the Milky Way and LMC, accreting at a rate above 
a few per cent of Eddington, show spectra peaking at 0.4-1.5 keV. 
They closely f ollow the standard L oc Tnmx dependence (except 
LMC X-1, see lOierliiiski & Dondliool and Fig.[3ji). The spread 
around it can result from varying contribution of the non-thermal 
emission. A shift of the microquasar G RO J 165 5 -40 to the ri ght 
from the L-T relation for a 7 M© BH JOrosz & BailviJll997t) is 
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probably related to a high spin of the BH there, which results in a 
higher disc temperature compared to the Schwarzschild BH. 

Two super-Eddington accretors in our galaxy, a persistent 
source SS 4 33 (.Dolan et al., 1 997) and a super-critical transient 
V4641 Sgr jRevnivtsev et alj 2002), show optical spectra with 
characteristic temperatures of 30 000-50 000 K, which we asso- 
ciate with the emission from the outer photosphere (and not with 
the s pherization temperature Tsp as proposed by Begelman et alJ 
l200d) . We estimate the accretion rate to be about lO'^ higher than 
the Eddington one. SS 433 is a bright UV sou rce with total luminos- 
ity of about 10*° erg s"* jPolan et aljl997l) . but is underluminous 
in the X-rays. V4641 Sgr, on the other hand, exceeded 12 Crab in 
the X-rays during the outburst and and was extremely bright in the 
optical band. Such a difference may be caused by a different mode 
of accretion and/or different inclination of the two systems. 

In Fig.[3jb) we present the data for ULXs with the luminosities 
exceeding lO'^^erg s^^ and the most reliable spectra obtained by 
XMM-Newton. The temperatures obtained bv lStobbart et al.l bood) 
from the fits with the blackbody and DISKBB models and the corre- 
sponding observed bolometric luminosities are plotted. A theoret- 
ical spectrum from a super-critical disc with a T{R) oc R~^^^ 
dependence, can, in principle, be represented as a sum of such 
components. We associate a high temperature component with the 
hot inner disc of Tmax~l keV. A high BH spin and an over- 
hea ting of the disc above the effective temperature at mo~2- 
20 ( lBeloborodovlll998l : ISuleimanov et alj|2002l. lKawaguchill2003h 
may be responsible for sometimes observed higher (up to 4 keV) 
temperatures. A soft, ~0.2 keV component may correspond to 
the spherization temperature implying the accretion rate rho — 
m-i/2(i.5j^/j;^p[keV])2«30-40 onto a stellar mass, 10-20 
Mq, BH. The observed higher luminosities can result from the ge- 
ometrical beaming. 

LMC X-1 deserves a more detailed discussion. Its apparent 
luminosity is sub-Eddington for a 10 M0 BH, but its track on the 
L-T diagram (see Fig. [3ji) is perpendicular to that for other BHs. 
This can be understoo d if its mass is only ^LSMp (i.e. it can be 
a neutron star, but see iHutchings et ai]|l983l:lEbrsawa et al.lll99lh 
and the source accretes at mo~3. The observed rather soft spec- 
trum is associated with the spherization region and results from 
obscuration of the central hot disc. In that case, however, we se- 
riously underestimate the true luminosity of the object. Assuming 
that it is five time larger than the apparent one, the behaviour of 
LMC X-1 on the L-T plane (Fig. [3)3) becomes consistent with that 
of the spherization temperature at rho ~l-3 for a 10 M0 BH. 



6 SUMMARY 

A black hole accreting at a super-Eddington rate is likely to produce 
strong winds (and jets) as observed in SS 433, the only known BH 
in the Milky Way accreting persistently at such a high rate. The 
bolometric luminosity can exceed the Eddington limit by a loga- 
rithmic factor ~ In mo, while the apparent luminosity can be factor 
of 5 higher because of the geometrical beaming. The observational 
appearance of such an object strongly depends on the inclination 
angle of the system to the line-of-sight. In edge-on systems, the 
central X-ray source is hidden by the wind and most of the radia- 
tion, with the characteristic temperature of ~ 10^-10^ K, escapes 
in the UV band. A face-on observer sees the hot inner flow with 
the flat in EFe spectrum extending from a few keV down to the 
temperature at spherization radius Tsp ~ rng^^'^ keV. 

Strong winds from the accretion disc explain naturally the 



presence around the ULXs of the expanding nebulae which are 
photoionized by the X-ray and UV radiation of the central source. 
An excellent agreement between the model and the data sup- 
port \2£wsJJFabrika&^Me^^ 2002; Fabrikal 
l2004lBegelman et alj2006l ; IVierdaYanti et al.l2006h that ULXs are 
super-critically accreting stellar-mass BHs similar to SS 433, but 
observed along the symmetry axis. 
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